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SUMMARY 

A lumped parameter model of the TF34 engine is formulated to study non- 
recoverable stall. Features of the model include forward and reverse flow, 
radial flow in the fan, and variable corrected speed. 

The purpose of the study is to point out those parameters to which 
recoverability is highly sensitive but are not well known. Experimental 
research may then be directed toward identification of the parameters in that 
category. 

Compressor performance in the positive flow region and radial flow in the 
fan are shown to be important but unknown parameters determining recover- 
ability. Other parameters such as compressor performance during reverse flow 
and in-stall efficiency have relatively small impact on recoverability. 

INTRODUCTION 

Compression system instability has received much attention over the past 
two decades. Intensive efforts in engine-inlet compatibility have done much 
to alleviate the problem of compressor stalls which was once quite serious. 

However, as the incidence of inlet-induced stalls was reduced by improved 
design, the tendency was then to operate the compressor at increased pressure 
ratio to achieve greater performance. The cost of this increased performance 
is a loss in stall margin which implies again increased stall incidence. As a 
result, the industry has taken the point of view that occasional compressor 
stalls are an acceptable price for improved performance. 

Many of the compressor stalls which occur recover, with only a momentary 
loss of power and without pilot action, and are of little consequence. Some 
stalls, however, cannot be cleared and are nonrecoverable. In these in- 
stances, the compressor settles into a stable operating point where it runs in 
rotating stall. Engine power output is very low, and internal temperatures 
are high enough to inflict damage. To recover from this condition, the engine 
must be shut down and restarted which requires deliberate pilot action and 
considerable time. Clearly, nonrecoverable stalls carry serious implications, 
and have been the object of intensive investigations in recent years. 

The predominance of compressor mapping effort has been devoted to the 
unstalled region. Literature reporting in-stall compressor performance is 
sparse (see, for example, references 1 and 2). Greitzer, et al., in ref- 
erences 3 to 5 present data for several machines and suggest modelling tech- 
niques to accommodate post-stall operation. 

The Air Force is currently contracting with General Electric and Pratt and 
Whitney to advance the technology in the nonrecoverable stall area. As a por- 
tion of their contract, GE has proposed to analyze TF34 data to be generated 
by LeRC. As a corollary to this TF34 experimental program, a post-stall model 
using characteristics specific to the TF34 has been developed. In this 




report, the model is described, and results in terms of factors pertinent to 
recoverability from stall are given. 

Description of Model 

A lumped parameter representation of the TF34 engine system was developed 
on a hybrid computer. Performance parameters are included so that post-stall, 
as well as normal behavior is simulated. To facilitate observation of the 
rapidly occurring events which follow stall, the model is operated 1000 times 
slower than real time. A block diagram of the engine system model is shown in 
figure 1. In the sections which follow, descriptions will be given of the 
manner in which each component of the engine system is represented. 

The symbols used in this report are defined in appendix A. A summary of 
equations is given in appendix B. 

Compressor 

The compressor is represented by overall maps, supplied by the manu- 
facturer, with pressure ratio a function of corrected airflow. Fluid inert- 
ance is lumped, and compressor volume is divided between upstream and down- 
stream stations. Mechanical speed is constant, but corrected speed varies. 

In the performance map given, compressor pressure ratio varies approx- 
imately as PNC^*^, (symbols are defined in appendix A) and corrected flow 
varies as PNC^. This is contrary to incompressible theory in which cor- 
rected flow varies directly with speed, and pressure ratio varies as speed 
squared. A possible explanation for these unusual exponents is that the 
compressor variable geometry is scheduled as a function of corrected speed in 
this steady-state map. 

The transients which are the subject of this study are rapid. It is 
judged that: 1) corrected speed cannot be sensed and computed that rapidly by 
the engine control system, and 2) the variable geometry actuators would be 
unable to respond quickly. Hence, the variable geometry would remain rela- 
tively fixed, and the reported map would not be valid. Accordingly, the speed 
line for PNC = 1. is used as a baseline. Other speed lines are obtained by 
shifting the PNC = 1 line by conventional incompressible theory; pressure 
ratio is shifted as PNC^, and corrected flow as PNC. 

In this approach, an in-stall pressure parameter \i>, proportional to 
downstream static pressure minus upstream total pressure divided by inlet 
density times wheel speed squared, is calculated as a function of a flow 
parameter proportional to flow velocity divided by wheel speed. Greitzer 
suggests a numerical value for the pressure parameter at zero flow based on 
observations of several machines. For non-zero flows, both positive and nega- 
tive, a parabolic relationship is assumed where ^ varied as and the 
proportionality constant is treated parametrical ly. 

The pressure parameter used is based on compressor downstream static pres- 
sure. Downstream total pressure is calculated by multiplying the downstream 
static pressure by a constant. Since compressor discharge mach number is 
always low, errors introduced by this simplification are minor. 

Hysteresis exists in the compressor stal 1 /recovery process. The operating 
point and the stalling flow are known. The value of flow at which recovery 
occurs is chosen to be greater than operating point flow, and is varied para- 
metrically. The compressor map which results is shown in figure 2. In this 
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figure, speed lines for PNC = 0.8, 0.9 and 1. indicate how performance shifts 
with corrected speed. 

When flow through the compressor reverses, the compressor inlet reference 
station is switched from the volume at the compressor inlet to the burner vol- 
ume. Compressor flow is then corrected by conditions in the burner. The cor- 
recting factor for compressor speed is not switched to burner temperature 
because that would result in an undesireable discontinuity at zero flow, and 
because data indicate that rotational speed is unimportant during reverse 
flow. 

Temperature performance in the in-stall, positive flow region is cal- 
culated from the pressure performance described above and hypothesized effi- 
ciency. Recent unpublished compressor rig data indicate a smooth transition 
in temperature ratio when going from unstalled to stalled operation. The 
stalled efficiency is based on the value required to achieve this smooth tran- 
sition in temperature ratio when going from unstalled to stalled operation and 
some low value at zero flow, with linear interpolation between. The zero flow 
efficiency is treated parametrically. Resulting temperature ratio 
characteristics, for PNC = 1. and several zero flow efficiencies are given in 
figure 3. 


Fan 

The fan is treated in a similar, but simpler, manner as the compressor. 
Unstalled performance is obtained from reference 6. Stalled performance is 
handled as in the compressor. For the fan, temperature ratio is assumed 
constant. 

Two fan maps, one for the hub and one for the tip are used, and are 
presented in figure 4. Radial crossflow between hub and tip regions is 
allowed, proportional to the pressure difference. Since the tip pressure 
ratio is slightly higher than the hub pressure ratio, a constant pressure 
source (linear map) is added in the hub to tip flow path. This is adjusted 
for zero quiescent cross flow. The flow impedance between the hub and tip 
regions is varied parametrically. 

Combustor 

The combustor is modelled simply. Heat is released as soon as fuel enters 
the combustor. During periods of reduced or reversed airflow, the mixture in 
the combustor becomes very rich because of the relatively slow fuel system 
dynamics. Instantaneous fuel-air ratio of the gases in the combustor is 
computed. If this ratio exceeds the stoichiometric ratio, only that portion 
of fuel for which the necessary air is present is allowed to react. 

Volume Dynamics 

Pressures and temperatures within the lumped volumes are computed using 
well established techniques (e.g., reference 7). For temperature determin- 
ation during reverse flow, the temperature multiplier is switched from the 
normal entering value to that of the volume(e.g., eq. B3). This has the 
effect of bleeding energy from the downstream volume and heating the inlet 
volume. 
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Rotational Speed 


The model does not include turbines or rotor dynamics; mechanical speed is 
constant. However, the reversed flow of hot gases raised the inlet temper- 
ature to the compressor and fan, which causes their corrected speeds to de- 
crease. This phenomenon is modeled. 

Control 

The fuel control is simply a constant Wf/P 3 , with no dynamics. Fuel 
pulse disturbances are introduced by programming a false P 3 , which is 
common test-cell practice. 

Compressor variable geometry is not modeled, as discussed earlier. 


Analytical Investigation 

The model is exercised to assess the importance of various parameters on 
the recoverability of the system following a compressor stall. Two procedures 
are required; one to stall the compressor and one to rate the stability of the 
system. 

Compressor stall is effected by fuel pulses. A false burner pressure sig- 
nal is applied to the simplified control wherein Wf/P 3 is constant. Fuel 
flow is ramped up to a value just sufficient to cause stall and is then ramped 
down to a minimum value. Equal times (0.010 sec.) are used for the ramp up 
and the ramp down. 

Turbine area is used as a stability rating criterion. At the time of 
compressor stall, turbine area is reset by some increment which can be posi- 
tive or negative. This increment is adjusted until the compressor just 
recovers. For the nominal case, this area is defined as S’. Then some param- 
eter is changed and a new recovery area identified by trial and error. A 
change in area is then defined as aA = Apg,^^. If the new area is larger 
than 7T, the parameter change has had a detrimental effect on recovery. 

Hence, positive aA/A is bad. 

Examples of transients which did and did not recover, plotted on a com- 
pressor map, are presented in figure 5. 

Results and Discussion 

A set of system parameters is chosen for a nominal case. The nominal 
case has positive recovery margin so that the turbine area required for 
recovery is less than its nominal value at the operating point. Many of the 
system parameters have to do with engine geometry and can be measured or cal- 
culated with reasonable assurance. Several parameters, however, are unknown 
and were necessarily chosen in an arbitrary manner. It is a purpose of this 
investigation to study the sensitivity of stall recovery to system parameters, 
especially those which are nebulous. In this way, experimental effort can be 
directed toward the evaulation of those parameters which are not well known 
but are nevertheless important. 

Some of these parameters will be discussed individually as warranted; the 
recovery sensitivities of all the parameters will be summarized. 
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Compressor In-Stall Performance 


In the in-stall positive flow region the pressure coefficient was modeled 
as i|»Q = KPC + 1.54 where is the flow coefficient, 1.54 is the value 

of i|)Q at zero flow suggested by Greitzer, and KPC is a parametrically 
varied constant. A similar relationship is used in the negative flow region, 
4«c = KNC <t>l + 1.54. 

Recovery is very sensitive to KPC (figure 6 ). The reason for this can be 
understood by referring back to figure 5. Results indicate that, for recov- 
ery, the trajectory should not intersect the in-stall characteristic in the 
positive flow region. The trajectory should remain below the in-stall char- 
acteristic at least until the flow at which recovery occurs is reached. (It 
must be kept in mind that PNC along the trajectory, and is generally not on 
one of the speed lines shown. For the recovering transient, PNC is approxi- 
mately 0.87 at recovery.) Since the slopes of the trajectory and the char- 
acteristic are similar an increase in the value of the characteristic would 
result in a significant decrease in the throttle area necessary to avoid 
i ntersection. 

Recovery is much less sensitive to the characteristic in the negative flow 
region, determined by KNC. The recovery sensitivities can be expressed as the 
slopes of the sensitivity plots in figure 6 at the origin. The values 

calculated from figure 6 are — a(aA/ A) _ — _ 35 ^ — a(aA/ A) — _ q 

a(aKPC/KPC) a(aKNC/KNC) 

A qualitative idea of this great difference in sensitivity is given in 
figjjre 7. In the positive flow region, the nominal characteristic is the 
aA/A = 0 curve. Above and below th^ nominal, characteristics corresponding 
to -10 percent and +10 percent aA/A, respectively, are shown. In the nega- 
tive flow region, the nominal characteristic is shown together with that for 
aA/A = +10 percent. It is clear that a minute change in the positive-flow 
stalled characteristic is equivalent to a large change in the negative-flow 
stalled characteristic. 


Fan Radial Flow 

Communication between the hub and tip regions of the fan has received much 
interest. Proximate splitters have been proposed and tried to prevent aug- 
mentor ignition pulses from entering the core. 

Flow in the radial direction is assumed to vary linearly with the pressure 
difference. A nominal value of flow resistance is arbitrarily chosen to be 
(0.01 Ibf sec/lbm in^). Since so little was known of the numerical value of 
this resistance, it is varied over a wide range. 

The recovery area required for a range of radial flow resistances is plot- 
ted in figure 8 . The curve is neither monotonic nor linear, but the general 
trend is toward less likely recovery as radial flow resistance is increased. 

Compressor Efficienty at Zero Flow 

Compressor efficiency is handled as previously described; arbitrarily 
chosen at zero flow and as required for a smooth temperature transition at 
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stall. The zero flow efficiency is found to have a relatively minor effect on 
recovery area. 


Fan Inlet Volume 

The forward boundary conditions are established by the volume at the fan 
entrance, and by the source impedance through which this volume communicates 
with the ambient. Since these parameters are unknown, they are varied over 
wide ranges. Recovery area is slightly affected by fan inlet geometry. 

Other Effects 

Table I summarizes the sensitivity of recovery area to all of the param- 
eters which were investigated. The results are tabulated as the slope of the 
function at the nominal point. The results are listed in order of decreasing 
sensitivity. A positive sign indicates that an increase in the value of the 
parameter has a detrimental effect on recovery, i.e., more turbine area is 
required for recovery. 

Again, one of the intents of this study is to point out areas which 
require further definition. In table I, some of the parameters high on the 
list, e.g., compressor discharge volume and inertance, involve engine geometry 
and may be readily evaluated. Others, e.g., compressor pressure performance 
at negative flow and compressor stal 1 /recovery hysteresis, are not well known 
but are relatively unimportant to the recovery process. 

Of concern are those parameters which are not known but are highly sensi- 
tive. Prime examples of these are compressor performance in the positive flow 
region and the resistance to flow between the fan hub and tip regions. 

Although these parameters will be difficult to measure experimentally, they 
are requisite for a meaningful engine system model. 

SUMMARY 

A lumped parameter model of the TF34 engine is formulated. Features of 
the model include forward and reverse flow through the machine, radial flow in 
the fan and variable corrected speeds. In-stall compression system perform- 
ance is necessarily hypothesized. Both recoverable and non-recoverable stall 
behavior are exhibited by the model. 

A study is conducted in which several system parameters are varied to 
determine their effect on recoverability. System stability is rated on the 
turbine area necessary to recover; a more stable system recovers with a 
smaller turbine area. 

The purpose of the study is to point out those parameters to which 
recoverability is highly sensitive but are not well known. Experimental 
research may then be directed toward identification of the parameters in that 
category. 

Examples of parameters sensitive but unknown are compressor performance in 
the positive flow region and radial flow in the fan. Other parameters such as 
compressor performance during reverse flow and in-stall efficiency have rela- 
tively small impact on recoverability. 
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APPENDIX A 
SYMBOLS 


area 

constant in linearization of Cp 
constant in linearization of Cp 
specific heat at constant pressure 
time increment 

performance function (i = 1,4) 

function 

enthalpy 

heat transfer coefficient 
heating value of fuel 
constant 

constant describing in-stall com- 
pressor performance at negative 
flow 

constant describing in-stall fan 
performance at negative flow 
constant describing in-stall com- 
pressor performance at positive 
flow 

constant describing in-stall fan 
performance at positive flow 
fluid inertance, length/(area g) 
rotational speed 
static pressure 
total pressure 

fractional corrected compressor 
speed 

fractional corrected fan speed 
heat flux 
gas constant 

flow resistance, a (pressure drop)/ 
A (flow) 
temperature 
time 
volume 
weight 
weight flow 
constant 

ratio of specific heats 
ratio of total pressure to standard 
atmospheric pressure 
incremental change 
efficiency 

ratio of total temperature to 
standard atmosphereic temperature 
flow parameter 
pressure parameter 


in^ 

”r2 


Btu/lbm 

J/kg k2 

Btu/lbm 

°R 

0/kg K 

Btu/lbm 

”R 

J/kg K 

sec 


sec 

Btu/lbm 


J/kg 

Btu/in^ 

sec °R 

J/m2 sec K 

Btu/lbm 


J/kg 

sec2/in2 


sec 2/ m2 

rev/min 

Ibf/in^ 


rev/min 

N/m2 




Btu/sec 

J/sec 

in Ibf/lbm “R 

N m/kg K 

Ibf sec/lbm in2 

N sec/ kg m' 

“R 

K 

sec 

sec 

in2 

m2 

Ibm 

kg 

Ibm/sec 

kg/sec 

- 

- 


Subscripts : 


0 

ambient 

1 

fan inlet 

2 

compressor 

2H 

fan hub discharge 

2T 

fan tip discharge 

3 

compressor discharge 

air 

air 

as 

air, stored 

ave 

average 

C 

compressor 

D 

design 

duct 

duct 

exit 

exit 

F 

fan 

f 

fuel 

fb 

burned fuel 

fs 

unburned fuel stored in combustor 

i 

initial 

H 

fan hub 

metal 

compressor metal 

S 

standard 

T 

fan tip 

X 

fan radial flow 


Superscripts; 

' denotes parameter switched during flow reversal 

Nominal 

* denotes quasi-steady parameter at compression component discharge 


Constants : 


KNC 

50. 

KNF 

5. 

KPC 

5. 

KPF 

0.5 

Rq 

1.644x10-3 

Rduct 

0.02715 

Rx 

0.01 


>• nominal values 
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APPENDIX B 
SUMMARY OF EQUATIONS 

NOTE: The terms within braces { } in the equations were switched when 

flow reversed. If flow was positive, in the normal sense, the upper term 
within the braces was used. The equations in which this applies are B3, Bll 
B20, B25, B31, B32, B35, B39, B40, B47, and B51. 

Fan Inlet 
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PNF = 


(B9) 




= 3.594x10 


-3 

PNF 


= { 


KPF 1 
KNF 


+ 0.11 


Fan Discharge 


Unstalled 


Stalled 


PNF “ 1\ PNF 


P2t/Pi - 


P* - 1 073 P I 1 + ’*'fi 


PNF ~ '’2\ PNF 


p* _ 1 073 P I 1 + '*'f 

^2H ■ *^1 ^ 0.9916 i 


T^-p = 1.11 


^2H “ ^1 


*^1 " (^2T ■ ^2jj *^1,1 


" Lu 


J ('^2H " ^2h) '^H,i 


(BIO) 

(BID 

(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(B17) 
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Fan Duct Inlet 



F~T " ^exit^ 
duct 


(“t - “x - "duct) * "2T,i 


-“xl iZ 1 - "duct <"T)2T,1 


2H 


^(WT) 

*2T " 


2T 


T = 
'2T 


M 


2T 


(B18) 

(B19) 

(B20) 

(B21) 

(B22) 


(WT)2h 


Compressor Inlet 



(B23) 

(B24) 

dt + (WT)2 hj 
(B25) 

(B26) 
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(WT)2h 

2H- 

(B27) 


(B28) 

^2 " ^ Jq ^^3 " ^3^ *^2,i 

(B29) 

^2H 

®2H ~ T. 

(B30) 


s 



(B31) 


(B32) 


PNC = 


^\V®2h), 


(B33) 


. o -,^-2 '^2V^/^2H 

( 1 )^ = 3.779x10 

(B34) 

4.. = 4 + 1.54 

(B35) 


KNC 
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Compressor Discharge 

Unstalled Stalled 


PNC^ ^ 3( PRC 


P3 = 1.026 P^j^ (1 + 0.8528 PNC'^ 

(B36) 


^3^^2H ^ ^ p / '^2V^^^2H 

~PN? PNC 


To = Tom I 1 + PNC^ — 


★ 

’3 ■ '2H 


(Pt/P,u)‘'' ■ - 1' 


n = + n 


^c = o 


(B37) 

(B38) 


Combustor 


as 


/ [{:i-n=w(«. -{!,))] 




f s 





dt + W 


fs, 1 


(B40) 


W 


G = 1 - 


f/a = 


a/f 

T5“ 


fs 


as 


(G > 0) 


(B41) 

(B42) 





(B43) 


Wfb = (1 - G) Wf 


(B44) 
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Q = 0.95 HV 


= 


«T3(«as * “fs> 


(MH). 


t * 

= V j - CpWjTj * Q j dt - (WH)j_, 


Ho = 


(WH), 


3 ■ Was W^s 


H. 


^3 - Cp 


^3 ^3 


ave 



(B45) 

(B46) 

(B47) 


(B48; 

(B49) 

(B50) 

(B5i) 

(B52) 
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TABLE I. - SENSITIVITY OF SMALL RECOVERY 
AREA TO VARIOUS SYSTEM PARAMETERS 


Parameter X 


Compressor discharge volume 

KPC, compressor positive flow pressure performance 

Compressor inertance 

Fan exhaust duct aP 

Fan tip discharge volume 

Fan inertance 

Fan radial flow resistance 
Fan hub discharge volume 
Fan inlet volume 
Compressor n at zero flow 

KNC, compressor negative flow pressure performance 
Fan inlet duct aP 
Compressor hysteresis 


a(aA/A) 

a(aX/X) 

^ = 0 
X 


- 1.67 

- 1.35 

+ 0.97 

- 0.83 

+ 0.56 

+ 0.38 

+ 0.24 

+ 0.20 

- 0.12 

- 0.09 

+ 0.05 

- 0.04 

+ 0.03 
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Inlet 

Ipss, 

»0 

Inlet 

volume, 

pi-fi 

Fan tip 
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map, Pjf Tjr 

Fan tip 
momentum, 

Wt 

Fan tip 
volume. 

Fan duct 

Iqss, 

'^DUCT 



Crossflow. 



Fan hub 
quasisteady 

Fan hub 
momentum. 

Fan hub 
volume. 

Compressor 
quasi steady 
map, P3, T3 

Compressor 

momentum, 

W2 

Burner 

volume, 

'’ 3,^3 

Turbine 

(choked) 

^3 



Figure 1. Block diagram of engine system model. 
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Figure 3. - Temperature performance maps for compressor with various zero- 
flow efficiencies. 



Figured - Pressure performance map. 
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(b) Fan tip. 

Figure 4. - Concluded. 



Figure 5. - Trajectories of recoverable and non recoverable transients. 
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Figure 6. - Variation in recovery area as a function of compressor in- 
stall performance coefficients. 
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Figure 7. - In-stall compressor map showing the effect of varying KPC and KNC. 




Figure 8. - Variation in recovery area as a function of radial flow 
resistance. 


